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We investigate the effects of a nearly uniform Bose-Einstein condensate (BEC) on the proper- 
ties of immersed trapped impurity atoms. Using a weak-coupling expansion in the BEC-impurity 
interaction strength, we derive a model describing polarons, i.e., impurities dressed by a coherent 
state of Bogoliubov phonons, and apply it to ultracold bosonic atoms in an optical lattice. We show 
that, with increasing BEC temperature, the transport properties of the impurities change from co- 
herent to diffusive. Furthermore, stable polaron clusters are formed via a phonon-mediated off-site 
attraction. 
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PACS numbers: 03.75.-b, 03.67.-a, 71.38.Mx, 71.38.Ht 



The lack of lattice phonons is a distinguishing fea- 
ture of optical lattices, i.e., conservative optical poten- 
tials formed by counterpropagating laser beams, and con- 
tributes to the excellent coherence properties of atoms 
trapped in them However, some of the most inter- 
esting phenomena in condensed matter physics involve 
phonons, and thus it is also desirable to introduce them 
in a controlled way into optical lattices. Recently, it has 
been shown that immersing an optical lattice into a Bosc- 
Einstcin condensate (BEC) leads to intcrband phonons, 
which can be used to load and cool atoms to extremely 
low temperatures Q • Here, we instead concentrate on the 
dynamics within the lowest Bloch band of an immersed 
lattice, and show how intraband phonons lead to the for- 
mation of polarons 0, Q]. This has a profound effect 
on lattice transport properties, inducing a crossover from 
coherent to incoherent hopping as the BEC temperature 
increases. Furthermore, polarons aggregate on adjacent 
lattice sites into stable clusters, which are not prone to 
loss from inelastic collisions. Since these phenomena are 
relevant to the physics of conduction in solids, introduc- 
ing phonons into an optical lattice system may lead to a 
better understanding of high-temperature superconduc- 
tivity [J [BJ and charge transport in organic molecules 
Q. Additionally, this setup may allow the investigation 
of the dynamics of classically indistinguishable particles 
0- 

Experimental progress in trapping and cooling atoms 
has recently made a large class of interacting many-body 
quantum systems accessible. For instance, the forma- 
tion of repulsively bound atom pairs on a single site has 
been demonstrated and strongly correlated mixtures 
of degenerate quantum gases have been realized [h| . In 
such Bose-Fermi mixtures, rich phase diagrams can be 
expected, including charge and spin density wave phases 
[ill . [T3 | , pairing of fermions with bosons .13j , and a su- 
persolid phase [l4|. Here we instead consider one atomic 
species, denoted as the impurities, confined to a trap- 
ping potential, for example an optical lattice, immersed 
in a nearly uniform BEC, as shown in Fig. [T] Based on 
a weak-coupling expansion in the BEC-impurity interac- 
tion strength, we derive a model in terms of polarons, 
which arc composed of impurity atoms dressed by a co- 
herent state of Bogoliubov phonons 0, 0] ■ The model also 
includes attractive impurity-impurity interactions medi- 
ated by the phonons 15,16[. An essential requirement for 




FIG. 1: (Color online) A quantum degenerate gas confined 
to an optical lattice is immersed in a much larger BEC. For 
increasing BEC temperature T, a crossover from coherent to 
diffusive hopping, characterized by J and E a , respectively, can 
be observed. The phonon-induced interaction potential Vi,j 
leads to the formation of off-site polaron clusters, separated 
by a gap Eb from the continuum of unbound states. 



our model is that neither interactions with impurities nor 
the trapping potential confining the impurities impairs 
the ability of the surrounding gas to sustain phononlike 
excitations. The first condition limits the number of im- 
purity atoms [Icj ]. whereas the latter requirement can be 
met by using a species-specific optical lattice potential 
[l7j |. Moreover, unlike in the case of self-localized impu- 
rities we assume that the one- particle states of the 
impurities are not modified by the BEC, which can be 
achieved by sufficiently tight impurity trapping. 

Model. — The Hamiltonian of the system is composed 
of three parts, H = H x + Hb + Si, where S x governs 
the dynamics of the impurity atoms, which can be either 
bosonic of fcrmionic. The BEC Hamiltonian Hb and the 
density-density interaction Hamiltonian Hi are 



Hb= J dr^(r) 



■^ + y e xt(r) + p t (r)^(r) 



Si = K J dr X t ( r )x( r )0 t ( r )0( r ) 5 

where x(r) is the impurity field operator and </>(r) is the 
condensate atom field operator satisfying the commuta- 
tion relations [^(r), <^(r')] = <5(r — r') and [</>(r), </>(r')] = 
0. The coupling constants g > and n account for the 
boson-boson and impurity-boson interaction respectively, 
mti is the mass of a condensate atom and I4 x t(r) a weak 
external trapping potential. Without yet specifying H x 
we expand x(r) = X^ 7 ^( r )^i where r] v (r) are a set of 
orthogonal mode functions of the impurities and a v (aj,) 
the corresponding annihilation (creation) operators, la- 
beled by the quantum numbers v. 
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A common approach to find the elementary excitations 
of the BEC in the presence of impurities is to solve the 
Gross-Pitaevskii equation (GPE) [19( for the full system, 
i.e., for k / 0, and to subsequently quantize small os- 
cillations around the classical ground state. To obtain a 
quantum description of the impurity dynamics, we in- 
stead solve the GPE without taking Hj into account, 
and express the BEC deformations around the impuri- 
ties as coherent states of Bogoliubov phonons. Specif- 
ically, we write ftr) = 4>o( r ) + 5(f>(r), with 0o( r ) = 
<f>Q (r) the solution of the GPE for k = 0. Provided 
that the impurity-boson coupling is sufficiently weak, 
i.e., \n\/gn (r)£, D (r) < 1, with £(r) = h/ y/m b gn (r) 
the healing length, n-o(r) — 0o( r ) an( i E> the num- 
ber of spatial dimensions, we expect that the devia- 
tion of (f>(r) from <f>o(r) is of order k, i.e., (Sftr)) oc k, 
where ( • ) stands for the expectation value. We in- 
sert <?!>o( r ) + 5(j)(r) into the Hamiltonian Hb + Hi, keep 
terms up to second order in n, and obtain the linear 
term k J dr x ( r )x{ r )^o( r )[^ft (r) + 6<f>(r)], in addition to 
the standard constant and quadratic terms in 5<f)(r) and 
5ft(r), since 4>o(r) is no longer the ground state of the 
system. 

In order to diagonalize the quadratic terms in 5<j) (r) 
and 5ft (r), we use the expansion Sftr) = 5Z^[u M (r)& M — 
v *i( r )^li]j where u M (r) and v^(r) are the solutions of the 
Bogoliubov-de Gennes equations [T3| for k — 0, and 
(bV) are the bosonic Bogoliubov annihilation (creation) 
operators, labeled by the quantum numbers \i. We as- 
sume that the mode functions rj v (r) are localized on a 
length scale much smaller than is set by 14 x t(r), and that 
/ dr 1 77^ (r) 1 2 1 ?7t- (r) | 2 w for v ^ r, i.e., the probability 
densities 1 77^ (r ) | 2 for different mode functions deviate ap- 
preciably from zero only within mutually exclusive spatial 
regions. In this case J dr </>o (r) [m^ (r) — u M (r)]^(r)^*(r) « 
and J drno(r)r] u (r)i]*(r) w hold for v ^ r, and 
hence the nondiagonal impurity-phonon coupling is neg- 
ligible. The total Hamiltonian can thus be rewritten in 
the form of a Hubbard-Holstein model [2(| H = H x + 

with huj^ the energies of the Bogoliubov excitations, the 
number operator h v = a\,a v , the dimensionless matrix el- 
ements M Vill = (n/hUfj,) J dr0 o (r) [M r ) ~ u /*( r )] l"^( r )| 2 
and the mean field shift E v = k J dr no(r)|77i,(r)| 2 . Wc 
obtain an effective Hamiltonian H e g including correc- 
tions to (j>o(x) of order k by applying the unitary Lang- 
Firsov transformation 0, 9 H eS = UHU\ with U = 
cxp [E„,„( M J,/# ~ M^bJhJ], which yields 

H cff = UHjj^ +^{E v -E v )n v -^E v h v {n v -l) 
^ V v , r n v h T + J~] fot^&t S M . (1) 

The transformed impurity Hamiltonian UH X U* is ob- 
tained using the relation UalW = cv v X\, where 
Xt = expfX^M*^}, - Af„ lM S M )] is a Glauber dis- 
placement operator that creates a coherent phonon cloud, 
i.e., a BEC deformation, around the impurity. In 



the limit where the BEC adjusts instantaneously to 
the impurity configuration, polarons created by a^Xi 
are the appropriate quasiparticles, and H c g describes 
a nonretarded interaction with the potential V UiT = 
hu,, (Af V)/1 M* „ + M* !fl M T)fl ) . The polaronic ' level 
shift E v = X M hu>^\M u ^\ 2 is equal to the characteristic 
potential energy of an impurity in the deformed BEC. 

We now turn to the specific case of bosons loaded 
into an optical lattice immersed in a homogeneous BEC 
0]. In the tight-binding approximation, the impurity 
dynamics is well described by the Bose-Hubbard model 

H x = -JT,(ij) aJaj+5^Ej«j(«3— 1 )+M < E J ;»j > whcrc 
fi describes the energy offset, U the on-site interaction 
strength, and J the hopping matrix element between ad- 
jacent sites [1, [2l[. The modes of the lattice atoms are 
Wannicr functions rjj (r) of the lowest Bloch band local- 
ized at site j, and denotes the sum over nearest 
neighbors. Noting that [£7, = 0, we find 

UH X W = - J E (XiarfXjaj + ^E%(% - 1) + » 
(hj) 1 3 3 

(2) 

with the corresponding matrix elements Afj iq = 
K\/noEq/ (Hujq) 3 /^(q), where q is the phonon momen- 
tum, e q = (h~q) 2 /2mb the free particle energy, ftw q = 
\/ e q( £ q + 2g n o) the Bogoliubov dispersion relation, and 
fifa) = ^ / dr|77j(r)| 2 exp(iq • r), with f2 the quan- 
tization volume. We note that, for |q| <C l/£, we have 
Mj q oc fj (q) / -\/|qf , whereas for |q| S> 1/^, one obtains 
M^ q «/ 3 -(q)/q a . 

The Hamiltonian H c g describes the dynamics of hop- 
ping polarons according to an extended Hubbard model 
[S|, provided that c ^> aJ/h, with c ~ \/gno/rnb the 
phonon velocity and a the lattice spacing. We gain qual- 
itative insight into the dependence of Vij and the con- 
stant polaronic level shift E v = E p on the system pa- 
rameters by considering a one-dimensional quasi-BEC in 
the thermodynamic limit. Wc assume a sufficiently deep 
lattice to approximate the Wannier functions by Gaus- 
sians of width a < f, and find V itj = {n 2 /^g) 
and E p = k 2 /2£g. We note that the interaction between 
impurities is always attractive. More importantly, for re- 
alistic experimental parameters, £ ~ a, and hence the 
off-site terms Vjj+i oc c~ 2a ^ are non-negligible. This 
interaction potential is a direct consequence of the lo- 
cal deformation of the BEC around each impurity, as 
shown in Fig. [T] For a set of static impurities at posi- 
tions Xj = aj, the overall deformation of the BEC den- 
sity to order k is given by n(x) = uq + £V (Xjb^bqXj) = 
"o-fc/ffOEje- 31 *"*' 1 ^- 

Coherent and diffusive transport. — We first consider 
coherent hopping of polarons at small BEC temperatures 
ksT -C E p , where incoherent phonon scattering is highly 
suppressed. Provided that £ = J/E p <g; 1, we can ap- 
ply the so-called strong-coupling theory [j| , and treat the 
hopping term in Eq. @ as a perturbation. Including 
terms of first order in £, we obtain the impurity Hamil- 
tonian 
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j TtiTlj , 



(3) 



with |Et = n+Kn -E p , U = U-2E p , and J = J ((XjXj)}. 
Here (( ■ )) denotes the average over the thermal phonon 
distribution and i,j are nearest neighbors. We find 
((XjX 3 )} = exp{-E q#0 |A/o, q | 2 [l-cos(q-a)](2iV q + l)}, 
with a the position vector connecting two nearest neigh- 
bor sites and Nq = (Q hul i/ k BT _ Thus, the hopping 
bandwidth decreases exponentially with increasing cou- 
pling constant k and temperature T. 

At high temperatures E p < k B T < k B T c (with T c the 
critical temperature of the BEC) inelastic scattering, in 
which phonons are emitted and absorbed, becomes domi- 
nant, and thus the transport of atoms through the lattice 
changes from being purely coherent to incoherent. We in- 
vestigate this crossover by deriving a generalized master 
equation (GME) for a single particle starting from H c g in 
Eq Q. Using the Nakaj ima-Zwanzig projection method 
[22|, we find that the occupation probabilities Pj(t) at 
site I and time t evolve according to the GME [23[ 



dPjjt) 
dt 



= / ds E i w u ^ (* - s ) - w u ^ (* - s )] • 

(4) 

where the effect of the phonon bath is encoded in the 
memory functions Wij(s), which are symmetric in (i,j). 
To first order in £, thus keeping only nearest-neighbor 
correlations, we find 



1(^1 



cxp{2^|A/ , q | 2 [l-cos(q-a)] 



x[(iV q + l)(e iw - s - l)+7V q (e- iw ^ - 1)]} 



(5) 



and Wij(s) = if i and j are beyond nearest neigh- 
bors. The nontrivial part of Wij(s) takes the values 
2(J/h) 2 at s = and 2(J/h) 2 in the limit s -» oo. In 
the regime k B T <C P p , we have J ~ J, and the mem- 
ory function Wij (s) is well approximated by 2 (J/h) 2 6 (s) 
[with 0(s) the Heaviside step function], which describes 
purely coherent hopping in agreement with For 
E p <C k B T -c fc s T c , we observe that J -C J, and Wij(s) 
drops off sufficiently fast for the Markov approximation 
to be valid, as illustrated in Fig. [^a). In this case one 
can replace p(i — s) by Pi(t) in Eq. (0J and after in- 
tegration over s the GME reduces to a standard Pauli 
master equation d t Pi(t) = ^ Wij[Pj(t) — Pi(t)] describ- 
ing purely incoherent hopping. The hopping rate is of 
the form w itj ~ J 2 exp(— E a / k B T) / (Hy/ k B TE a ) [HH, 
where E a ~ E p is the activation energy. 

The temperature-dependent crossover from coherent to 
diffusive hopping in a quasi-one-dimcnsional (ID) system 
is apparent in the evolution for a time r ~ 1/wij of 
a particle initially localized at lattice site j = 0. The 
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FIG. 2: (a) Memory function h 2 W itj {t) /{2J 2 ) versus time for 
k B T = (dotted line), 2 E p (dashed line) and 10 E p (full line). 
It drops off rapidly for k B T 3> E p , indicating the dominance 
of incoherent hopping, (b) Coefficients hA/J (dashed line) 
and hy/B/J (full line) versus temperature, obtained from the 
numerical solution of the GME in Eq. {J4j for an evolution 
time r = Wh/J « 0.8 x 10~ 3 s. The condition A = \fB is 
satisfied at k B T/E p « 2.3, with E p /k B « llnK. The lattice 
(wavelength A = 790nm) contains a single 41 K atom with 
J = 2.45 x 10~ 2 E R , k/E r \ = 2.3 x 10~ 2 , the recoil energy 
E R = (2-rrh) 2 /2m a \ 2 , and m a the mass of the lattice atom. 
The BEC consists of 87 Rb atoms with no = 5 x 10 6 m _1 and 
g/E R \ = 8.9 x 10" 3 . 



mean-squared displacement of the lattice atom, l 2 (t) — 
l 2 Pi(t), can be decomposed as l 2 (t) = At+Bt 2 into in- 
coherent and coherent contributions, characterized by the 
coefficients A and B, respectively. The crossover takes 
place when incoherent and coherent contributions to I 2 (r) 
are comparable, i.e., \[B = A. Fig. [2b) shows A and \/B 
as functions of T, where l 2 (t) was obtained by numeri- 
cally solving the GME using the memory function Wij 
in Eq. ©. We find that, for a 41 K- 87 Rb system [H] un- 
der standard experimental conditions, the crossover takes 
place well below the critical temperature of the BEC (see 
caption of Fig. [2]). 

Polaron Clusters. — We now discuss the formation of 
polaron clusters for k B T < E p , based on Hamiltonian 
H {1) in Eq. ©, and assume that the bosonic impuri- 
ties are in thermal equilibrium with the BEC. At these 
temperatures Vij, and J are well approximated by their 
T = values. We consider the limit U ^> Vjj+i, U ^> J, 
and adiabatically eliminate configurations with multiply 
occupied sites. Keeping only nearest neighbor interac- 
tions, we obtain approximate expressions for the bind- 
ing energy Eb(s) m (s — l)V, J+ i of a cluster of s po- 
larons located in adjacent sites and the lowest energy 
band E k (s) w -E b {s) - 2J s (V h j+i) 1 ' 8 cos(fca) [H, with 
k the quasimomentum. This band approximation is in 
good agreement with the results from exact diagonaliza- 
tion of iP 1 ) using the full interaction potential Vij, as 
shown for three polarons in Fig. [3Ja). 

This model predicts a decreasing average cluster size 
with increasing temperature. For a small system with 
N = 3 polarons, we calculate the probability of finding 
a three-polaron cluster P3 = *^2 v (n u ), where the sum is 
taken over all states with energies e v < E g {2) + E g (l), 
Eg (N) is the ground state energy of an A^-polaron cluster, 
and the occupation probabilities are given by the Boltz- 
mann law (n v ) oc exp(— e v /k B T). Analogously, we deter- 
mine the probability P2 of finding a two-polaron cluster 
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FIG. 3: (Color online) (a) Energy spectrum of three polarons 
in a ID optical lattice. The solid line shows the lowest energy 
band Es(k) characterizing off-site three-polaron clusters. The 
lattice (wavelength A = 790nm, M — 31 sites, periodic bound- 
ary conditions) contains 133 Cs atoms with J = 7.5 x W~ 3 Er, 
U = 50E P , and k/E r \ = 1.05 x 1CT 1 . The BEC consists 
of S7 Rb atoms with n = 5xl0 6 m _1 and g/E R \ = 4.5xl0~ 2 . 
(b) Probability of finding a three-polaron cluster (solid line) or 
a two-polaron cluster (dashed line) versus temperature. The 
parameters are q/ErX = 6.5 x 10~ 2 , k/Er\ = 1.32 x 10 _1 , 
U = 2.2E P , M = 27, and the rest as for (a). 



Conclusion. — We have demonstrated that the dynam- 
ics of bosonic impurities immersed in a BEC is accurately 
described in terms of polarons. We found that spatial co- 
herence is destroyed in hopping processes at large BEC 
temperatures while the main effect of the BEC at low 
temperatures is to reduce the coherent hopping rate. Fur- 
thermore the phonons induce off-site interactions which 
lead to the formation of stable clusters which are not af- 
fected by loss due to inelastic collisions. Using the tech- 
niques introduced in this paper qualitatively similar phe- 
nomena can also be shown to occur for fermionic impu- 
rities. In cither case these effects can be controlled by 
external parameters and lie within the reach of current 
experimental techniques. 
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FIG. 4: (Color online) (a) Density-density correlation in a sys- 
tem of three polarons indicating the formation of off-site three- 
polaron clusters, (b) Corresponding momentum distribution 
(n k ). k/ErX = {4.0, 6.1, 8.1, 10.1, 12.1} x 10 -2 where a 
higher value corresponds to a more localized correlation in 
(a) and a more spread distribution in (b). The other parame- 
ters are g/E R X = 6.5 x 10~ 2 , M = 27, U = 3E P , and the rest 
as in Fig. [3ja). 

or bipolaron. The results are shown in Fig. [Hb). The 
probability of having a three-polaron cluster goes down 
with increasing temperature and for T = E p /ks ~ 18nK 
is essentially zero for the parameters chosen. For this 
three-polaron cluster, three-particle loss is negligible due 
to the on-site repulsion U. Decreasing the value of U 
gives a significantly increased P2 compared to the values 
shown in Fig. |3fb). 

The clustering of polarons leads to their mutual expo- 
nential localization. This is illustrated by the density- 
density correlations (fiihi+j) for the three-polaron clus- 
ter in Fig. 0Ja). With increasing attractive interaction 
Vij, the mutual localization gets stronger, leading to an 
increased broadening of the momentum distribution, as 
shown in Fig. BJb). This allows polaron clusters to be 
identified in time of flight experiments. We note that the 
transition from a superfluid to a Mott insulator also leads 
to broadening of the momentum distribution as, e.g., ob- 
served in [l0[. However, using Bragg spectroscopy @ 
would allow the unambiguous distinction of boson clus- 
tering from this transition. 
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